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ABSTRACT 



Aims. We study horizontal supergranule-scale motions revealed by TRACE observation of the chromospheric emission, and investi- 
gate the coupling between the chromosphere and the underlying photosphere. 

Methods. A highly efficient feature-tracking technique called balltracking has been applied for the first time to the image sequences 
obtained by TRACE (Transition Region and Coronal Explorer) in the passband of white light and the three ultraviolet passbands 
centered at 1700 A, 1600 A, and 1550 A. The resulting velocity fields have been spatially smoothed and temporally averaged in order 
to reveal horizontal supergranule-scale motions that may exist at the emission heights of these passbands. 

Results. We find indeed a high correlation between the horizontal velocities derived in the white-light and ultraviolet passbands. The 
horizontal velocities derived from the chromospheric and photospheric emission are comparable in magnitude. 
Conclusions. The horizontal motions derived in the UV passbands might indicate the existence of a supergranule-scale magneto- 
convection in the chromosphere, which may shed new light on the study of mass and energy supply to the corona and solar wind at 
the height of the chromosphere. However, it is also possible that the apparent motions reflect the chromospheric brightness evolution 
as produced by acoustic shocks which might be modulated by the photospheric granular motions in their excitation process, or ad- 
verted partly by the supergranule-scale flow towards the network while propagating upward from the photosphere. To reach a firm 
conclusion, it is necessary to investigate the role of granular motions in the excitation of shocks through numerical modeling, and 
future high-cadence chromospheric magnetograms must be scrutinized. 
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1. Introduction 

Quasi-steady convective flows of different scales have been sug- 
gested to play an important role in the processes of mass sup- 
ply and energy transport across different layers of the solar at- 
mosp h ere (e.g.. lFoukall|197l iKrijger et alJl2002t iMarsch et all 



2008; iDammasch et alJ l2008tTCurdt et alJ I2008I) . As is well 
known, granulation with a typical 1-Mm size and 5 -minute life- 
time are ubiquitous in high-resolution photospheric i mages (e.g., 
iTitle et alJfl98l iBerrilli et al.ll200l I.Tin et alJl2009h . Groups of 
theses granules tend to move in a systematic way, which is char- 
acterize d by a cellular convec tive motion at scales of the order of 
32 Mm dLeighton et al.ll 962). This larger-scale flow was termed 
supergranulation, with cells of which the boundaries coincide 
with the chrom ospheric network and lanes of magnetic con- 
centrations (e.g. . [ Leigh ton et al][l~9 62: Simon & Leighton 1964; 
ISimon et al.lfl988HWang & Zirinll 19881) . 

Under the assumption that the granules can be consid- 
ered as tracers of t he underlying larger -scale velocity fields 
( Simo n et al. | [1981 iRieutord et all [2001) of magnetoconvec- 
tion, one can derive the related horizontal flow speed by 
the widely used LCT (Lo cal correlation trackin g) tec hnique 
dNovember & Simon|[l988l) . Recently. IPotts et all (120041) devel- 
oped a new method for tracking flow fields. This so called ball- 



Send offprint requests to: H. Tian 



tracking method is more noise-tolerant and highly efficient for 
analysing large data sets. The flow field inferred by this method 
has a similar accuracy to that obtained by LCT. 

The horizontal motions in the chromosphere have not been 
scrutinized and are not understood well. Few attempts have 
been made to measure chromospheric proper motions, e.g., 
by ap plying LCT to the Ha images obtained in active re- 
gions Wi & Molownv-Horasll995l:IChae et al.ll2000l;|Yang et alJ 
120031) . By using this method, the typical horizontal velocity in 
the chromospheric network was found to be 1000-1500 m/s 
(lYi & Molownv-Horasl ll995). However, the small-scale struc- 
tures present in Ha images appear to be elongated and thus 
are different from the cellular-shaped granules. Moreover, the 
Ha line is very sensitive to and strongly influenced by dynamic 
events. Thus, the explanation of the horizontal velocities as mea- 
sured from Ha images is not straightforward. 

Besides Ha, the Ca n H & K lines have also been exten- 
sively used to explore dynamics in the chromosphere. The chro- 
mospheric network, which coincides with the magnetic network 
and outlines the supergranular boundaries, is the most prominent 
structure on images of Ca n, while the cell interiors (internet- 
work regions) which are enclosed by the network occupy most 
of the quiet-Sun area. Internetwork areas are filled with inter- 
mittent grainy brightenings. These emission features are usually 
termed H(K)2,, cell grains because their intensities peak dramat- 
ically just blueward of the Ca n H & K line centers. These inter- 
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network grains (or cell grains) are also present on wider-band (a 
few A) Ca n filterg rams, but at a lower contrast and a slight phase 
shift dRutten et al.lll999allbh . The sizes of internetwork grains are 
about 1 ~4 Mm, and their lifetime s are usually less than 12 min- 
utes dRutten & Uitenbroek|[T99Tl: iTritschler et al]|2007l) . For a 
detailed review of the Ca n H2„ & K-2 V cell grains, we refer to 
iRutten & Uitenbroekl(ll99ll) . 

By comparing the quiet-Sun images o btained by TRACE 
(Transition Region and Coronal Explorer. lHandv et alJ fl999) 
in its three ultraviolet (UV) passbands (centered at 1700 A, 
1600 A, and 1550 A) and a co-temporal Ca n K filtergram, 
IRutten et alJ dl999al) found a high degree of spatial agreement 
between the Ca n K image and the TRACE images, in particular 
the 1700 A one. They concluded that the three TRACE chan- 
nels portray internetwork grain phenomena as well as Ca n K 
filtergrams do. Therefore, TRACE image sequences provide an 
excellent means to study the dynamics of internetwork grains 
(IRutten et al.lll999ah . 

By applying here the balltracking technique for the first time 
to the TRACE images obtained in the passband of white light 
and the three UV passbands centered at 1700 A, 1600 A, and 
1550 A, we have made an attempt to study horizontal photo- 
spheric and chromospheric motions in a large area of the quiet 
Sun. Our analysis revealed a striking correlation between the 
horizontal velocities derived in the white-light passband and the 
UV passbands, a finding which might shed new light on our 
understanding of mass supply to the corona and solar wind. 
Alternatively, our finding may provide insights into the excita- 
tion and propagation of the shock sequences which might pro- 
duce the bright internetwork grains. 

2. Data analysis and results 

The TRACE image sequence analysed here was obtained on 14 
October 1998 from 08:20 to 09:30 UTC in a quiet-Sun region 
around disk center. Images in the four passbands as mentioned 
above were taken with a cadence of 21 s. The sizes of the images 
are 512x512 pixels squared, with a pixel size of 0.5". 

The standard software for reducing TRACE data was applied 
to this data set, including the removal of cosmic rays, subtraction 
of dark currents, normalization of the counts, and so on. Then 
we used the IDL procedure RotJCY.pro available in SolarSoft to 
take into account the solar rotation effect. The satellite jitter was 
subsequently removed with an accuracy of 0.025", by applying 
the cross-correlation technique to enlarged (by a factor of 20) 
images. For computational efficiency, we resampled the data to 
obtain a cadence of 42 s, and then extracted a sub-region with a 
size of 200" x 200" from the coaligned images. 

Each image was spatially Fourier filtered to remove solar os- 
cillations and observational noise and thus to rem ove the large 
scale variations, as described in iPotts et a fl (12001 . In addition, 
the 1700 A, 1600 A, and 1550 A images have the problem that 
the intensity variation is very non-linear, with a small number of 
very bright points that tend to dominate the filtered data. These 
data were scaled by raising the intensity to a power of 0.4 to 
compress the intensity data into a more uniform range, before 
filtering. We kept the Fourier components with wavelength be- 
tween 1" (the minimum measurable wavelength) and 3.5" (up- 
per limit of the granule/grain size). Figure Q] presents the first 
images of the sequences taken with the passbands of white light 
and 1700 A. It is clear from the figure that the size of the chro- 
mospheric grains inside cell interiors is comparable to that of 
the photospheric grains. The brightness range of the cell grains 



seems to be larger than that of the photospheric grains. The most 
prominent difference between the photospheric emission and the 
chromospheric emission is the much more enhanced network 
emission in the chromosphere. But the network bright points 
only constitute a minor part of the area, and thus their contri- 
bution to the tracked velocities is minor. Figure Q] also shows the 
corresponding filtered images, which reveal clearly that both the 
dark and bright points inside the cell interiors are comparable in 

size (or width). 

In the balltracking technique as developed by IPotts et aTl 

(2004), one considers a photospheric image in a three- 
dimensional representation by regarding the brightness as cor- 
responding to a geometrical height. Small floating balls are 
dropped on to this surface and tend to settle in the local min- 
ima between some adjacent granules. As the granulation pattern 
evolves, the balls will be pushed around and thus follow and re- 
veal the motion of the granules. It has been demonstrated that the 
results from th is method have a n accuracy similar to those pro- 
duced by LCT dPotts et al.l2004l) . This method has been success- 
fully applied t o high-resolution continuum images obtained by 
MDI/ SOHO dPotts et al.|l2007l: IPotts & Diver! 120081: llnnes et all 
120091) and SOT/Hino de lUttie et alJl2009h . 

As mentioned bv lPotts et alJ d2004l) . the balltracking method 
is not limited to flows in the photosphere, but will also work for 
any velocity field in which there are visible moving features of 
known scale length. Here we made the first attempt to apply this 
method to TRACE images in the passband of white light and 
several UV passband s. We set the radiu s of the tracking ball at 
1.2". As described in lAttie etail d2009t) . the output of the ball- 
tracking procedure are velocity fields derived from the displace- 
ments of many individual tracking balls, which reflect the fast 
and stochastic small-scale granular motions. To extract the un- 
derlying large-scale velocity field th e data needs to be spa tially 
smoothed and temporally averaged. iRieutord et alJ d2001l) sug- 
gested that the velocity field may not be faithfully described by 
granular motions at spatial scales less than 2.5 Mm, or at tem- 
poral scales shorter than half an hour. Our tracking results have 
been smoothed over 3.5 Mm and running-averaged over the en- 
tire 70 minutes. The tracking routine actually tracked the data 
twice, once with the normal data and once with the inverted data, 
and gave very similar results in both cases. We took the average 
of the two results. The error estimate for the tracking was ob- 
tained from the difference of these two results. The root mean 
square values of the errors are well below those of the tracked 
velocities in all of the four passbands. 

Finally, since the velocity values obtained by balltracking 
are usually under estimated as compared to the true velocities 
dPotts et al] I2004T) . we need to divide the balltracked veloci- 
ties by a scaling factor. Starting with the de-rotated and jitter- 
removed data cube, we chose a sub-area from the first frame of 
the data, and offset each subsequent frame in turn by a small dis- 
tance, representing an additional imposed velocity. Four veloci- 
ties were chosen to produce four offset data cubes, which were 
then tracked using balltracking. By applying a linear fit to the 
mean velocity derived from each offset data cube versus the off- 
set velocity, we obtained a scaling factor in each passband. The 
derived velocities were divided by this value. Figure|2]shows the 
result of the fitting in each passband. The values of the slopes 
are used as velocity calibration factors. 

The distributions of the resulting velocities are presented in 
Fig. [3] It turns out that the velocity values are comparable in 
different passbands. The median values of the derived veloci- 
ties are 247 m/s, 258 m/s, 258 m/s, and 284 m/s, in the pass- 
bands of white light, 1700 A, 1600 A, and 1550 A, respectively. 
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White light-before filtered 




White light-after filtered 
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Fig. 1. Images taken at 08:20 in the passbands of white light (upper) and 1700 A (lower). The original and filtered images are 
presented on the left and right, respectively. The size of each image is 200" x 200". 



Following the method in lPotts et al.l d2004i) . we estimated an un- 
certainty of 29 m/s for the velocity measurement. 

Figure [4] presents the velocity fields for the four passbands. 
The cell boundaries, which were calculated by ap plying the auto- 
matic recognition algorithm developed by iPotts & Di ver (2008) 
to velocities obtained in the corresponding passband, are shown 
as the background (lanes). A direct comparison of the cell 
boundaries as calculated from the velocity field in white light 
and from that in the passband of 1700 A with the 1700 A inten- 
sity image was made and is shown in Fig. [5] 

We also calculated the linear Pearson correlation coefficients 
of each component of the velocity field as well as the veloc- 



ity magnitude, between each pair of passbands. The results are 
listed in Table Q] The average correlation coefficients of the un- 
filtered intensity images between each pair of passbands are also 
listed therein. 

3. Supergranular horizontal flows in the 
photosphere 

Photospheric horizontal flows have been extensively investigated 
through ground based observations. Image sequences obtained 
by the TRACE satellite are suitable for motion tracking in the 
photosphere as they offer long uninterrupted runs over large ar- 
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Fig. 2. Calibration factors obtained by a linear fit to the mean 
velocity derived from each offset data cube versus the offset ve- 
locity, in the passbands of white light (black), 1700 A (red), 
1600 A (green), and 1550 A (blue). The values of the slopes 
are used as velocity calibration factors. 
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Fig. 3. Histogram of the measured horizontal velocities in the 
passbands of white light (black), 1700 A (red), 1600 A (green), 
and 1550 A (blue). 

Table 1. The linear Pearson correlation coefficients of the inten- 
sity (/), two components of the velocity field (v x and v y ), as well 
as the velocity magnitude (v) between each pair of passbands. 



passbands 


/ 


v x 


Vy 


V 


WL/1700 


0.13 


0.63 


0.71 


0.35 


WL/1600 


0.09 


0.58 


0.62 


0.26 


WL/1550 


0.06 


0.50 


0.61 


0.22 


1700/1600 


0.96 


0.81 


0.81 


0.63 


1700/1550 


0.82 


0.73 


0.71 


0.48 


1550/1600 


0.92 


0.78 


0.81 


0.61 



eas (up to 512" x 512") with no seeing problems. The spatial 
resolution of 1" is just high enough to resolve the granular struc- 
ture. 

By applying the balltracking technique to white light images 
observed by TRACE, we have reconstructed the well-known 
photospheric supergranular flow pattern in a quiet-Sun region 
with a size of 200" x 200". From Fig. [5] we find that the calcu- 
lated cell boundaries generally match the network pattern as seen 
in the 1700 A image. In the 1700 A image, there are still a few 
segments of the bright lanes which are not reproduced from these 
calculations. It might be due to the effect of temporal variation, 
height variation, UV contamination of the white-light emission, 
smoothing and averaging of the tracking results, or the fact that 
small granules below the detection limit of TRACE can not be 
traced properly. 

The v elocity magnitude i n w hite light is consis tent with the 
results in Rrijg er et al.l(l2002|) an d lAttie et alj d2009l), but s maller 
than those in lWang et all (Il995l) and lRoudier et all <ll999h . This 
should be the result of the different smoothi ng and averaging of 
the tracked velocity fields. As mentioned in Atti e et al.l d2009l) . a 
weaker smoothing and averaging will increase the velocity val- 
ues. Our choice of the smoothing and averaging scales is a com- 
promise between a more accurate velocity measurement and a 
high er resolution to dis tinguish flows with sharp gradient. 

iKriigeret ail (120021) applied the LCT technique to a TRACE 
white-light image sequence in a small region, and concluded that 
it is possible to use the TRACE observations to measure photo- 
spheric horizontal velocity fields. Here we have demonstrated 
that the supergranular flow pattern in the photosphere can be re- 
covered by applying the highly efficient balltracking method to 
the white-light images observed by TRACE in a large area of the 
quiet Sun. 

4. Horizontal motions revealed through the 
chromospheric emission 

As mentioned previously, the horizontal flow field in the 
chromosphere has not been understood well. Only a few at- 
tempts have been made to measure chromospheric proper mo- 
tions by applying LCT to Ha images obtained in active re- 
gions dYi & Molownv-Horasll99"5l:IChae et al.ll2000l:lYang et all 
2003). The TRACE observations allow a direct comparison of 
possible large-scale horizontal motions in the chromosphere 
with those in the underlying photosphere. 

From Fig.|4]we can see that the velocity fields obtained in the 
three UV passbands also show cell structures which are similar 
to those in the white-light passband. The diverging flow-like pat- 
terns in some cells are highly coincident in different passbands. 
The correlation coefficients listed in TableQ]are large, again indi- 
cating a high correlation between the derived horizontal motions 
in different passbands. The cell boundaries calculated from the 
velocity field in the 1700 A passband also coincide more or less 
with the chromospheric network, as shown in Fig. [5] According 
to our knowledge, this is the first time that supergranulation- 
scale horizontal motions, both in and above the photosphere, are 
revealed in a large area of the quiet Sun. 



4.1. Formation height of emission from TRACE UV filters 

According to lHandv etall (1 19991) and lWorden et al.l d 19991) . the 
emission in the 1700 A passband (4-10xl0 3 K) is UV contin- 
uum originating around the temperature minimum; the emis- 
sion in the 1600 A passband (4-10xl0 3 K) has contributions 
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Fig. 4. Velocity vector fields as calculated by using the image sequence in the passband of white light (upper left), 1700 A (upper 
right), 1600 A (lower left), and 1550 A (lower right). The lanes represent cell boundaries as calculated from the velocity field in the 
corresponding passband. The size of each image is 200" x 200". 



from the UV continuum and several lines like Fe n and C i, but 
the dominant contributor is the UV continuum from the region 
of temperature minimum. Accordi ng to the VAL-C solar atmo- 
sphere model dVernazza et alJ[l98Th . the temperature minimum 
region is located at about 500 km above Tsm= 1 , corresponding to 



the upper photosphere/lower chromosphere i 


Worden et al.ll999: 


Rutten & Kriiaer 2003: McAteer et alJl2004 


). 



The emission from the 1550 A passband is mainly a com- 
binat ion of C iv and the underlying UV continuum emis- 

119991). As mentioned in 



sion dHandv et al J 1 19991 : IWorden et al 



iKriigeret alJ d200ll) andlRutten et al 



(1999a), the C iv lines 



dominate the emission of the 1550 A passband in active regions, 
while in the quiet Sun contributions from the UV continuum and 
several weak lines (e.g. C i lines) dominate in the 1550 A pass- 
band. jHandxetal] d 1999b mentioned that the 1550 A image re- 
sembles an image taken in the center of the Ca n K line with a 
bandpass of a few Angstroms. So the major contributor to the 
emission of the 1550 A passband seems to be the UV contin- 
uum, and thus the 1550 A passband should be dominated by the 
chromospheric emission. 
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Fig. 5. The yellow lanes representing cell boundaries as calculated from the velocity field in white light (left) and from that in 
the passband of 1700 A (right) are superposed on the averaged intensity image in the 1700 A passband. The size of the image is 
200" x 200". 



Similar to many papers such as iRrii ger et al.l d2001l) and 
iMcAteer et ail (T2004). we use the term "chromosphere" even for 
the upper photosphere since images of the TRACE UV pass- 
bands largely differ from the photospheric images and essen- 
tially sample the chromospheric network. We realize that there 
are limitations to assigning the height of formation to emis- 
sion from the broad band-pass filters on the TRACE UV chan- 
nels. However, by comparing the quiet-Sun images obtained by 
TRACE in three UV passbands (centered at 1700 A, 1600 A 
and 15 50 A) and a co-temporal Ca n K filtergram, iRutten et al.l 
(1999a) concluded that there is a high degree of spatial corre- 
spondence between the Ca n K image and the TRACE images, 
and that the three TRACE channels reveal internetwork emission 
features as well as Ca n K filtergrams do. Similar conclusion was 
also reached by lHandv et al.l d!999l) . Thus, we believe that the 
structures visible on the TRACE UV filters are well formed at a 
certain height, probably the lower part, of the chromosphere. 



The cross-talk (overlap) between emission from different 
passbands might also influence our results. First, since our data 
were obtained shortly after TRACE was launched, the transmis- 
sion c urves of UV passbands presented in Fig. 9 of lHandv et al.l 
(1999) should have not changed too much, and thus the leak- 
ing of the white-light signal into the UV emission should be 
marginal. Second, there is a large part of overlap between emis- 
sion of the 1700 A and 1600 A passbands, which accounts for the 
very large correlation between velocities in the two passbands. 
This overlap does not change our conclusion, since our main 
finding is the correlation between velocities in the white-light 
and UV passbands, not the correlation between two UV pass- 
bands. Third, the UV contamination of the white-light emission 
might not be trivial and is likely to be one of the reasons why 
some bright chromospheric network lanes are not reconstructed 
by the calculated cell boundaries. This is, however, counteracted 



by a higher correlation between the velocities than the intensi- 
ties. 



4.2. Nature of internetwork grains 

Internetwork grains (or cell grains) have been extensively stud- 
ied through ground bas ed observations using C a n H & K lines 
(for a review, see Rutt en & Uitenbroekl 1199 ll) . As mentioned 
above, IRutten et al.l d!999ab found a high degree of spatial cor- 
respondence between the Ca n K image and the TRACE UV 
images (passbands of 1700 A, 1600 A, and 1550 A). They 
concluded that the three TRACE channels portray internetwork 
grain phenomena as well as Ca n K filtergrams do. 

Internetwork regions occupy most of the quiet Sun area. So 
with the balltracking technique we are mainly tracking the mo- 
tions of the internetwork grains in TRACE UV filters. Before in- 
terpreting the observed horizontal motions in the chromosphere, 
we need to understand the nature of these internetwork grains. 

In fact, there are long-lasting disputes on this issue. Direct 
evidence for the spatial correspondence b etween the cell grains 
and internetwork mag netic elements (iLivingston & Harvevl 
1 19751) was provided by ISivaraman & Livingston! ([1982). The 
ma gnetic origi n of the cell grains was further confirmed 



Hi gra 



Dame (1985), Dame & Martic (198"; 



Sivaraman (1991), and Sivaraman et al. 



(2000). Based on Hinode observations, Ide Wiin et alj (120081) 
found that some of the magnetic elements are associated with the 
chromospheric cell grains. Recently, although no obvious corre- 
lation between the magnetic flux density and the Ca n H inter- 
network brightness was obtained, lYang et ail J2009) were still 
inclined to believe that there should be a correlation between the 
two, and they attributed the lack of correlation to the low cadence 
of Ca ii H images they used. 
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However, others doubt this conclusion. And through ob- 
servations and modelings they believe that the occurrence of 
internetwork grains does not depend on magnetism but is 



through the low chromosphere ( Rutten & Uitenbroek 


1991 


Remling et al. 1996; Carlsson & Stein 


1997; Lites et al. 


1999 


Loukitcheva et al.l2009t).IWorden et al.l( 


1999) compared the cell 



grains observed in TRACE UV passbands and the corresponding 
photospheric magnetic field, and suggested that the internetwork 
magnetic field is essentially uninvolved with the production of 
cell grains in the 1600 A passband, but could have a small in- 
volvement with the cell grains in the 1216 A (Lya) passband. 

iBrandt et all d 19921 11994 found that there are two types of 
cell grains: (i) oscillatory K2,. flashers, appear only a few times 
with a modulation of 3-minute oscillation, brief spatial mem- 
ory, 5-10 grains per cell; (ii) persistent K-line grains, retain- 
ing long-term identity while traveling through the cell, flash- 
ing with 3-5 minutes periodicity, about 1 grain per cell. The 
short-living grains were suggested to be generated by acoustic 
shocks and have nothing to do with the magnetic field. On the 
contrary, a spatial correspondence between internetwork mag- 
netic features and persistent flashers was observed. The per- 
sistent flashers probably represent the chromospheric signa- 
tures of newly-emerged strong-field ephemeral regions on their 
ways to be part of the network. This find ing was confirmed 
bv lNindos& Zirinl dl998h . iLites et alJ dl999h also found persis- 
tent flashers but pointed out that they are v ery rare. Whi l e with 
UV image sequences obtained by TRACE, Rrii ger et al J (12 00 1 ) 
identified many such long-living features, and McAteer et al. 
(2004) also found that ~50% cell grains are persistent flash- 
ers j^ersj^teji^flashers were called internetwork bright points 
bv Ide Wiin et al.l d2005), in which they were suggested to cor- 
respond to the strongest internetwork magnetic elements. 

The question of whether or not the chromospheric cell grains 
are correlated with the magnetic flux in internetwork regions is 
very important to study the chromospheric heating and dynam- 
ics. Here we don't intend to answer whether or not the magnetic 
field plays a role, or a major role, in the production of chromo- 
spheric cell grains, since there are still debates on this issue and 
it is not so eas y to distinguish betw een the two possibilities. As 
pointed out by iKamio et al.l d2006l) . we need more detailed ob- 
servations to solve this problem. Instead, we discuss the impli- 
cations of our results and point out directions of future studies, 
with respect to each of the two suggested natures of cell grains. 

4.3. Magneto-convection in the chromosphere 

As shown in Fig. |4] the measured horizontal velocity fields in 
the chromosphere and photosphere show a similar supergranule- 
scale cell-like pattern, which prompts us to think if they repre- 
sent an upward extension of supergranular horizontal flows to 
the chromosphere. 

The idea that supergranular flows may penetrate into higher 
layers has already been suggested by several authors. It is known 
that motions of granules in the photosphere follow the super- 
granulation plasma flows because these plasma flows are ob- 
served in Doppler shift. In the middle chromosphere (forma- 
tion height of Si 11 A 18 17), such Doppl er shift was also investi - 
gated and demo n strated to be existing dNovember et al.|[l979t) . 
iNovember etaP (1 19791) even predicte d that supergranu l ar ve- 
locities should be evident in the TR. November et al. (119821) 
found that large-scale patterns of up- and downflows in the mid- 
dle chromosphere correlate well with those seen in the photo- 



sphere and concluded that vertical flows of supergranular scale 
appear to extend into the chromosphere. Height variation of 
the velocity fie ld in the phot osphere has also been studied. 
Observations by Deubnei] (Il97ll) indicated that the vertical mo- 
tions increase with height, and that the horizontal supergranular 
flow decreases slightly, in the photosphere from the formation 
height of C 1 ,15380 to that of Mg 1 /15173. The latter is formed 
at the tem perature minimum (upper photosphere/ l ower chro- 
mosph ere) dNovember et al.lfl982tlNovemberlll989l) . lGiov anelli 
dl980 ) found that the supergranular horizontal velocities show 
no significant variation with height over the range of formation 
of C 1 A9111, Fe 1 /18688, and Mg 1 A8806, but there is a sub- 
stantial reduction to abo ut one-ha l f of th is at the level of Ca 11 
/18542. More recently, iDel Morol d2007l) investigated the 3-D 
photospheric velocity field of a supergranular cell, and proposed 
that large flow features can penetrate into the upper photosphere. 

ISimon & Leighto n (1964) suggested that the supergranula- 
tion flow field carries the internetwork magnetic field to cell 
boundaries, forming the pattern of the chromospheric network. 
So far many observations have confirmed that most internetwork 
magnetic elem ents move toward the boundaries o f the magnetic 
netwo rk (e.g.. iMartinl 11988b ISimon et alJ 11988b IWang & Zirinl 
1988; Zhang et al. 1998). Using magn etograms obtained b y the 
Big Bear Solar Observatory (BBSO), IZhang et alJ d2006l) esti- 
mated the emergence rate of the new flux (ep hemeral reg i ons) to 
be ~5xl0~ 19 cirT 2 day 1 in the quiet Sun. IWang et all d 19951) 
calculated the mean area of internetwork magnetic elements, 
which is about 10" 2 . Then the mean size and height of an in- 
ternetwork magnetic element should be about 3". According to 
IVernazza et alJ d 198 lb . the temperature minimum, which marks 
the boundary between the photosphere and chromosphere, is 
located at 500 km above T5oo=l- So in the internetwork re- 
gions, there should be a large number of magnetic elements 
which extend above the photosphere. And when these internet- 
work magnetic elements are swept by the supergranular flows 
to the magnetic network, we should see the horizontal motions 
of their footpoints and apices in the photosphere and chromo- 
sphere, respectively. Unfortunately, high-quality magnetograms 
of the quiet Sun above the photosphere are still difficult to ob- 
tain. Thus, currently it is difficult to validate the existence of 
these strong internetwork magnetic el ements a b ove t he photo- 
sphere. However, in an active region iLin et a 11 (120061) found a 
correlation between localized emission enhancement above the 
photosphere (in the chromosphere and TR) and a moving mag- 
netic feature, which fo l lowed the boundary of a supergranulation 
cell. IZhang & Zhangl (UpOO) analyzed quiet-Sun photospheric 
and chromospheric magnetograms observed by the vector video 
magnetograph at Huairou Solar Observing Station, and found 
that all visible variations in the photosphere had corresponding 
variations in the chromosphere, although the H/3 line, which they 
used to derive the chromospheric magnetogram, are complicated 
and ma y include signals from both the photosphere and chromo- 
sphere. ISolanki et alJ d2003l) observed a set of rising magnetic 
loops and found that magnetic signatures of the loop legs are 
clearly present in the magnetogram of the u pper chromosphere. 
Recen tly, based on Hinode observations, Shi moio & Tsuneta] 
(2009) even found a clear correlation between coronal activities 
and photospheric minority magnetic polarities in the polar coro- 
nal hole. All these observations seem to suggest that horizon- 
tal movements of internetwork magnetic elements should reveal 
horizontal motions both in and above the photosphere. 

If the chromospheric internetwork grains have a 
one-to-one correspondence to enhanced magnetic fields 



dSivaraman & Livingstonl[T982l; iDame et al 1 ll984HDamelll985b 
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Dame & Martidll987l Il988t ISivaramad[l99ll: [Sivarama rTeTall 
2000), or the persistent flashers which are related to strong 



internetwork magnet i c elements ( Brandt et al .l 119921 1994 
Nindos & Zirinlll998b iKriiger et alJl200lt iMcAteer et al.l l2004: 
de Wiin et al J 12005) play a dominant role in the chromospheric 



internetwork emission, the velocity fields revealed by TRACE 
UV observations naturally represent the horizontal component 
of magneto-convection in the chromosphere. By comparing 
the velocity fields obtained in white light with those in the 
three UV passbands, we can conclude that the supergranular 
pattern does not completely disappear above the photosphere. 
The diverging flows in some cells are still very clearly seen in 
the chromosphere. Some supergranular motions are extending 
into the overlying formation layers of the 1700 A, 1600 A, 
and 1550 A emission. The comparable velocities derived in 
the photosphere and chromosphere are also easy to understand 
because they reflect horizontal motions of the same features, 
although at different heights. 

The existence of a supergranule-scale magneto-convection 
in the chromosphere may shed new light on the study of mass 
and energy supply to the corona and solar wind at the height 
of the chromosphere. It is believed that part of the magnetic 
network flux opens into the corona in the shape of funnels, 
whilst the rest of the network consists of a dense population 
of low-lying loo ps with lengths less than 10 4 km and vary- 
ing o rientations (IDowdv et all 1 19861: iMarsch & full 1997b iPeterl 
2001). Magnetic loops reaching to the height of the chromo- 
sphere and TR can be swept by the supergranular flow from 
the cell interior to its boundaries, where they can interact with 
magnetic funnels, and by reconnection may supply mass and 
energy to the funnels dAxford et al.l 11999b lAiouaz et alj [2 005; 
Aiou az|l2008t iHe et all 120071: iMdntosh et al.l 120071: iTian et al l 
2008a, 2009). Reconnection at the interface between cool side 
loops and the network flux tubes may occur in the chromosphere 
and TR, result ing in the solar wind outflo w dTu et al.l 12 005: 



He et al.H 2008: Buc hner & Nikuto wski 2005; Tian et al. 2008b, 
2010) or upflows along loop legs (Tia n et alj|2008al 12009). and 
downflows at lower layers. Thus, horizontal motions represent- 
ing supergranule-scale magneto-convection in the chromosphere 
and TR are vital in such processes as mass supply and energy de- 
livery to the corona and solar wind through funnels. Our finding 
may represent the observational evidence for these motions. 

However, we should point out that the extension of strong 
magnetic elements into the chromosphere can only be explic- 
itly confirmed when the magnetic field above the photosphere 
can be accurately measured, and that the existence of such 
supergranule-scale magneto-convection in the chromosphere 
can only be confirmed by following the evolution of the chromo- 
spheric magnetic field with high cadence and high resolution. 



4.4. Chromospheric brightness evolution induced by 
photospheric granular motions 

As mentioned above, some authors doubt the magnetic origin 
of the chromospheric cell grains, but believe that cell grains 
are caused by weak acoustic shoc ks which propagate upwards 
through the low chromosphere (|Rutten & Uitenbroekl 199l(: 
Reml ing etal]|1996l: ICarlsson & Steinlll997t ILoukitcheva et al l 
2009 ULites et al.l(ll999h confirmed the existence of magnetism- 
related grains but pointed out that they don't play a dominant 
role in the chromospheric internetwork emission. If most cell 
grains are indeed of acoustic origin, then the horizontal motions 



we derived by using TRACE UV observations are likely to repre- 
sent chromospheric brightness evolution induc ed by the shocks. 

In the one-dimensional simulation of ICarlsson & Steinl 
( 1997), a piston located at the bottom of the computational do- 
main (100 km below T5oo=l) drives waves through the atmo- 
sphere, then propagating waves near or just above the acoustic 
cutoff frequency interfering with higher frequency waves that in- 
duce them to steepen rapidly and form shocks near 1 Mm above 
T5oo=l- The chromospheric shocks produce a large source func- 
tion, yielding the high emissivity of internetwork grains. If this 
scenario is correct, it is essential to investigate the horizontal 
distribution of pistons at the bottom layer, which w as not ad- 
dressed in the simulation of ICarlsson & Steinl (1 19971) due to its 
one-dimensional restriction. 

Several att empts have been mad e to investigate this issue. As 
pointed out by Rutt en et al.l dl999bl) . turbulent convection seems 
to supply th e pistons that excite sh ock sequences and produce the 
cell grains. Rimmele et al. I (fl99l found that enhanced acoustic 
wave activities, the so-called "acoustic events", occur preferen- 
tially in dark intergranul ar lanes. Thi s spa tial corre spondence 
was later confirmed by | Goode et all d 19981) and IStrous et al.l 
d2000h . lHoekzema & Rutted d 19981) and lHoekzema et alJdl998l) 
further pointed out that dark intergranular lanes tend to show ex- 
cess ~3-minute waves in the photosphere. These transient acous- 
tic waves are suggested to be excited by small granule s that un- 
dergo a rapid collapse dRastl ll999: Ska rtlien et al.ll2000l) . 

The relationship between enhanced wave activities in the 
photosphere and enhanced chromospheric emission was also 
investigat ed by several autho rs. For example, through wavelet 
analyses, Kam io et alj (|2006) found that the occurrence of chro- 
mospheric brightenings is correlated with enhanced 5 mHz ve- 
locity oscillations in t he chromosphere and the photosphere. 
iHoekzema et alj ([1998) concluded that the preferential align- 
ment between ~3-minute waves and dark photospheric inter- 
granular lanes does not survive to the chromospheric heights; in- 
stead, there seems to be a correspondence between excess chro- 
mospheric brig htness and intergranula r lanes at a time delay 
of 2.5 minutes. Hoekzema et al T d2002l) found that sites of en- 
hanced wave activity in the granulation preferentially co-locate 
with exceptionally bright chromospheric internetwork grains, at 
a delay of about 2 minutes which might represent the sound 
travel time to t he chr omosphere. While through a statistic study, 
ICadavid et all d2003l) found that 72% of the G-band darkening 
events are followed by an enhanced chromospheric emission 
2 minu tes later; in the rema ining 28% cases, the timing is re- 
versed. ICadavid et all J2003) also found that the G-band darken- 
ing events are usually accompanied by transient enhancement of 
the measured magn etic field , indica ting collapse of intergranular 
lanes. According to iGoodel d2002l) . a collapse in the intergran- 
ular lanes can produce upward-propagating waves which subse- 
quently lead to chromospheric brightenings. 

The horizontal motions we derived in the chromosphere 
show a similar pattern to the photospheric supergranulation. If 
most chromospheric cell grains are produced by acoustic shocks, 
our results should indicate that the acoustic shocks are mod- 
ulated by photospheric or subsurface flows. Wedemever et al.l 
(2004) pointed out that convection motions play an important 
role in the excitation of acoustic waves. It is well known that 
photospheric granules tend to move in a systematic way as char- 
acterized by the supergranulation. If the collapse of intergranu- 
lar lanes is really the source of acoustic waves that produce the 
chromospheric cell grains, we may expect a coherence between 
the photospheric supergranulation and the motion of the collaps- 
ing site of intergranular lanes. In other words, the collapsing site 
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of intergranular lanes has a tendency to move with the granules. 
This will also naturally yield a comparable horizontal velocity 
in the photosphere and chromosphere. 

As pointed out by Hoekzem a & Ruttenl dl998l) . it is worth- 
while to use photospheric flow tracking to enable studies of 
chromospheric response to photospheric or subsurface excita- 
tion sites while following migrations of the latter over mesoscale 
and larger distances. Our study focuses on the comparison be- 
tween horizontal motions in the photosphere and chromosphere, 
and thus provide some insights into this issue. The coherent 
supergranule-scale behavior between motions in the two layers 
seems to indicate that the chromospheric cell grains mark loca- 
tions where acoustic events follow on granular collapses in the 
evolving intergranular lanes. 

It is also possible that the propagation of the shock waves is 
influenced by the supergranulation. The shock wave flux might 
be advected by the supergranular flows toward the network while 
it is propagating upward from the photosphere and dissipated in 
the chromosphere. 

However, the resolution of the TRACE data we used here 
is not very high. More high-resolution and high-cadence ob- 
servations of different layers are needed to further study the 
dynamics of granular ev olution in intergran ular lanes and the 
chromospheric response (iRutten et al.l l2008h . And detailed 3- 
dimensional numerical simulation should also be done to inves- 
tigate the role of granular motions in the excitation of acous- 
tic shocks and the subsequent production of chromospheric cell 
grains. 

5. Summary and conclusion 

We have applied the highly efficient balltracking technique to 
TRACE images obtained in the white-light band and three UV 
passbands centered at 1700 A, 1600 A, and 1550 A. We have 
demonstrated that the supergranular flow pattern in the photo- 
sphere can be recovered by applying this tracking method to the 
white-light images observed by TRACE. This is the first time 
that horizontal motions in a large quiet area of the solar chro- 
mosphere have been investigated. Our analysis revealed a strik- 
ing correlation between the horizontal velocities derived in the 
white-light band and the UV passbands. 

The interpretation of our finding is not straightforward, since 
we tracked the apparent motions of the chromospheric internet- 
work (cell) grains, the nature of which is still under debate. If 
the cell grains (or most of them) correspond to enhanced in- 
ternetwork magnetic elements, the velocity fields revealed by 
TRACE UV observations should represent the horizontal com- 
ponent of magneto-convection in the chromosphere. Then our 
finding seems to provide evidence for the way on which mass 
and energy are supplied to the corona and solar wind at the 
height of the chromosphere, which is predicted or suggested by 
many recent observational and modeling studies. However, as 
believed by many authors, the cell grains may be entirely pro- 
duced by acoustic shocks propagating upward to the chromo- 
sphere. If the cell grains are indeed of acoustic origin, the veloc- 
ity fields revealed by our TRACE UV observations should reflect 
the motion pattern of the short-living chromospheric brightness 
as induced by acoustic shocks. Then the striking correlation be- 
tween the horizontal velocities derived in the UV and white-light 
passbands seems to indicate that the excitation of enhanced wave 
activity is modulated by the photospheric granular motions, or 
that the shock waves are advected by the supergranular flows to- 
ward the network while it is propagating upward from the pho- 
tosphere and dissipated in the chromosphere. 



We conclude that it is important to investigate the role of 
granular motions in the excitation of shocks through numerical 
modeling. In addition, future high-resolution and high-cadence 
observations including dopplergrams, magnetograms, and imag- 
ing of the photosphere and chromosphere are needed to investi- 
gate and understand the coupling between the two layers. 

We realize that our results and conclusions need to 
be checked and confirmed by future observations and with 
improved techniques. The mission of the Solar Dynamic 
Observatory (SDO), which was launched in February 2010, will 
provide full-disk high-resolution photospheric magnetograms 
and chromospheric images. It is likely that these data will better 
our understanding of the evolution of the chromospheric emis- 
sion. 
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